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Abstract 

For spherically symmetric relativistic perfect fluid models, the well-known 
Buchdahl inequality provides the bound 2M/R ^ 8/9, where R denotes the 
surface radius and M the total mass of a solution. By assuming that the ratio 
p/ p be bounded, where p is the pressure, p the density of solutions, we prove 
a sharper inequality of the same type, which depends on the actual bound 
imposed on p/p. As a special case, when we assume the dominant energy 
condition p/p ^ 1, we obtain 2M/R ^ 6/7. 
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1 Introduction 



1 Introduction 

In General Relativity, the masses M of equilibrium configurations with a given 
fixed radius R are limited. Indeed, from the theory of black holes, it is clear that 
2M/ R < 1 is required, since otherwise apparent horizons must necessarily appear. 
In 1959, H.A. Buchdahl f¥j proved a stricter inequality: 

H^£<? (1) 

Buchdahl's inequality holds for spherical stellar models whose matter is ordinary 
perfect fluid matter satisfying a non-decreasing barotropic equation of state p = 
p{p) relating the (positive) energy density p and the (positive) pressure p. The 
major appeal of Buchdahl's inequality lies in the fact that it leaves a gap between 
the maximally attainable value of 2M/R and the threshold of horizon formation. 
Buchdahl's inequality is a classical result in General Relativity and a discussion of 
it is included in almost every textbook. Its significance in astrophysics is due to 
the fact that ([1]) yields a bound on the surface red shift (gravitational red shift) of 
stellar models. 

There exists several bounds on 2M/ R under different assumptions than Buchdahl's. 
In [5], the inequality ([1]) is shown to hold also for anisotropic matter, where the 
tangential pressure is assumed to be less than the radial pressure. In p], collisionless 
matter is considered instead of fluid matter; in addition, a very general inequality 
is derived that is independent of the specifics of the matter content. Recently, the 
authors of [7j considered a large variety of assumptions. In particular, they obtain 
inequalities by imposing merely the dominant energy condition; see also the remarks 
in the conclusions. 

In this paper we return to Buchdahl's original realm, i.e., we investigate ordinary 
perfect fluid matter that satisfies a non-decreasing equation of state. Although 
the Buchdahl inequality ([1]) is sharp under these assumptions, the solution that 
maximizes the quantity 2M/ R is rather unphysical: The stellar model satisfying 
2M/R = 8/9 consists of an incompressible fluid, i.e., p = const throughout the 
entire configuration, and p diverges as r — > 0. Accordingly, for this limiting solution, 
p/ p is unbounded, so that in particular the dominant energy condition is violated. 
In this paper we derive a generalization of Buchdahl's inequality that applies if — in 
addition to Buchdahl's original assumptions — we also assume a bound on p/p; this 
new inequality is formulated in Theorem 13. li Its major property can be condensed 
into the following key statement: The more restrictive the bound on p/ p, the more 
restrictive the bound on 2M/R. The extremal cases covered by Theorem 13.11 are 
the case p/ p ^ oo (i.e., no bound), where the original Buchdahl inequality ([1]) is 
recovered, and the case where p/p ^ e for an arbitrarily small e > (which is the 
most restrictive bound), which leads to 2M/R ^ const e. When we assume the 
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2 Dynamical systems formulation 



dominant energy condition (i.e., the bound p/ p < 1), we find that 

2M 6 

Since the methods of our proof are based on the dynamical systems approach to 
relativistic perfect fluid models, which is due to [6], we begin by briefly introducing 
these techniques. 



2 Dynamical systems formulation 



Static spherically symmetric relativistic perfect fluid solutions (relativistic stellar 
models) are solutions of the relativistic equations of hydrostatic equilibrium: 

47rr^/9 , (3a) 

\ pc^ J \ mc^ J \ &r J 

Equation (pb|) is the famous Tolman-Oppenheimer-Volkoff (TOV) equation [U [TO] . 
The function m = m{r) is the so-called mass function; when R is the surface radius 
of a stellar model, we denote by M = m(R) its total mass. A solution of ([3]) gives 
rise to a static spherically symmetric perfect fluid spacetime, which is represented 
by the metric 

where the "relativistic potential" v is given by du/dr = —{dp/dr)/{p + pc?). Note 
that r is the standard areal radial coordinate. To form a determined system, these 
equations are supplemented by a barotropic equation of state p{p). 

We assume that p and p are non-negative and related by a barotropic equation of 
state p{p) that is sufficiently smooth for p > Q. (In principle, equations of state 
that are only piecewise smooth do not pose any further difficulties; however, for 
simplicity, we confine ourselves to the simplest case.) In the following we choose 
units so that G = \ and c = 1. 

Although the system ^ is adequate for several purposes, it is not suitable for prob- 
ing generic features of the solution space associated with large classes of equations 
of state. For example, in the standard proof of Buchdahl's inequality, see 0] or [3], 
the system is replaced by an alternative system that is based on a choice of variables 
different from (m,p, r). To remedy the defects of the TOV equation, in [B], the sys- 
tem ([3]) is reformulated as a dynamical system. This reformulation turns out to be 



dm 
dr 
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dr 
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2 Dynamical systems formulation 



the key to the derivation of several quahtative results such as theorems on the finite- 
ness of R and M, and results on qualitative features of the M-i2-diagram. In the 
following we shall present a brief summary of this dynamical systems formulation, 
since it constitutes the basic ingredient in the proof of our main theorem. 

The two-dimensional system ([3|) allows a reformulation in terms of a regular au- 
tonomous three-dimensional system of differential equations defined on a compact 
state space; see [6]. To obtain this dynamical system formulation one makes a 
variable transformation from (m,p > 0, r > 0) to the variables 



m 




uj = uj{p) =p ; (5) 



in principle, uj{p) can be a more or less arbitrary function that is strictly mono- 
tonically increasing and sufficiently smooth. For our purposes, however, the choice 
UJ = p is sufficient. 

Since the variables are unbounded, we replace these variables by bounded 

variables {U,V,Q) G (0, 1)^ by simply defining 

l + U l + V 1+UJ ^ ' 

In addition, we introduce a new independent variable A, which is a measure of r, by 

J^ = {i-urHi-vr'. (7) 

alogr 

Using these definitions turns ([3]) into a dynamical system on the bounded state space 
X = (0, 1)3 = {{U, V, n)\U,V,n e (O, l)}, which is the unit cube. The dynamical 
system reads 

^ = U(1- urn -V){3- AU) -TH] (8a) 

^ = V{1- V)[{2U -l){l-V + 2aV) + {l-T) H] (8b) 
aX 

^ = -n{i-n)H, (8c) 

CIA 

where H = H{U,V,n) = {1 + a)V{l - U + aU) . (8d) 

The r.h. side contains two functions that encode the information on the equation of 
state: 

T = - — , o- = - . (9) 

p dp p 

In the context of the system ([8]), these functions are to be viewed as functions of 
O (instead of functions of p). For details on the derivation of the equations dSD we 
refer to l6l. 
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3 A GENERALIZATION OF THE BUCHDAHL INEQUALITY 



Given an equation of state, the TOV equation admits a one-parameter family of 
regular solutions, which are solutions of ([3]) with a regular center of spherical sym- 
metry. The parameter characterizing this family is most conveniently chosen to be 
the central pressure of the solution (i.e., pc = p|r=o)- (For a proof see 0; the dy- 
namical systems analysis of L6J provides a simple alternative.) Since m ~ ATrr^pc/S 
in the limit r — > (where pc = p{pc) is the central density), we find that regular 
solutions satisfy 

^7^^, V ^0, n^Qc (10) 

as r — > 0; in addition, 

X^-oo ^ r^O. (11) 

Therefore, in the dynamical systems picture, a regular perfect fluid solution is repre- 
sented by an orbit ([/(A), y(A),$7(A)) in the state space X that emerges from what 
we call the hue of regular initial values {{U, V,n)\U = 3/4, V = 0,n £ {0, 1)}. We 
conclude by noting that the limit A — > oo corresponds to the limit r ^ R, where R 
is the surface radius of the perfect fluid solution; this is because Q — > (and thus 
p ^ 0) is this limit; for details, see [6]. 

Remark. In the definition of the variable V, the factor 1 — 2m /r appears. This does 
not pose any problems, since this factor is always positive for relativistic perfect 
fluids; cf. p]. One can even show that it is not necessary to resort to additional 
considerations to ensure that 1 — 2m /r > 0; this issue can be dealt with completely 
within the dynamical systems formulation: If 1 — 2m /r > holds initially (as is true 
for regular solutions), then 1 — 2m/r > for all r corresponding to A G (— cxd,C)o) 
and thus in particular for all r < R. This is simply because 1 — 2m /r can be written 
as 1 — 2m /r = (1 — V){1 — V + 2aV), which is positive, since < V < 1 for all 
A S (—00,00). To show that also 1 — 2M/R > (which corresponds to the limit 
r ^ R, or, equivalently, to A — > 00), a straightforward dynamical systems analysis 
of the (j-limits of solutions in X is sufficient; see [6j. 



3 A generalization of the Buchdahl inequality 

Consider the function 

, _ , 9M \2{2M +V)- (37W + Vf] , , 

$ = (^(M V) = — — ^ — 12) 

^ (1-27W)(3A^ +P)2 ' 



where 

For each constant 



M = — , V = A-Kr^p . (13) 
r 



5G[0,4), (14) 
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3 A GENERALIZATION OF THE BUCHDAHL INEQUALITY 



the inequality 

(15) 

defines a Buchdahl type inequality. This is simply because it is an inequality relating 
m/r with the pressure p. 

Example. In the special case g = 0, the inequality (fTSj) reduces to 

2{2M+V)-{3M+Vf -;?0; (16) 
alternatively, when we use standard variables, we obtain 

^ ^ ^ (l - Gnr^p + Vl + QTTr'^p) , (17) 

which is the (complete) Buchdahl inequality. Consider a regular perfect fluid solution 
(m(r),p(r)) that satisfies this inequality; let R denote the surface radius, which is 
determined by the vanishing of the pressure, and M = m{R) the total mass of the 
solution. Setting P = in ([TB]) . or p = in (flTl) . which corresponds to evaluating 
the inequality at the surface i2, yields the so-called Buchdahl surface inequality 

2M 8 

The Buchdahl inequality is due to H. A. Buchdahl [4]; it holds for all regular perfect 
fluid solutions that obey an equation of state that is non-decreasing. 

Setting "P = in (fT5]) (for an arbitrary value of 5 E [0,4)) corresponds to evalu- 
ating this inequality at the surface of a perfect fluid solution. A simple algebraic 
manipulation shows that <I'|-p=o ^ 5 is equivalent to J^\r ^ (4 — 5)/(9 — 2g) or 

2M 8-2g 

^ . (19) 

R 9-2g ^ ^ 

The inequality (jlSp thus naturally incorporates a surface inequality of the type (jl9p . 

Theorem 3.1. Consider an equation of state p = p{p) that is non- decreasing; let 
a = pI p. Every associated regular perfect fluid solution that satisfies 

a<^-l (20) 
V9 

for some constant g € [0,4), obeys the inequality ^ ^ g and consequently satisfies 

2M 8-2g 

^ . (21) 

R 9-2g ^ ^ 
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4 Proof of the theorem 



Remark. In the special case g = 0, the assumption (j20p disappears. The theorem 
then reduces to the statement that every regular perfect fluid solution satisfies the 
Buchdahl inequality. The Buchdahl inequality is thus naturally included as a special 
case of the theorem; note in addition that the proof of the Buchdahl inequality 
becomes particularly simple in our set-up. 

Considering the special case g = 1 m the theorem leads to an interesting corollary: 

Corollary 3.2. Consider an equation of state p = p{p) that is non- decreasing. Every 
associated regular perfect fluid solution that satisfies the dominant energy condition 
(i.e., a = p/p < 1) obeys the inequality <I> ^ 1 and consequently satisfies 

2M 6 

In the subsequent section we give a proof of the theorem which is based on the 
dynamical systems formulation of the problem. 



4 Proof of the theorem 

We begin by proving that, in the dynamical systems description, regular solutions 
satisfy U ^ 3/4 (where U — > 3/4 as A — > — oo), see also [6]. This statement is almost 
trivial, see the subsequent remark; however, we give a rather illustrative proof which 
already makes use of those methods that will turn out to be pivotal in the proper 
proof of Theorem 13. 1[ 

Remark. In standard variables, the statement U ^3/4 reads p{r) ^ pir), where 
p = 3m/(47rr^) is the "average density". This inequality is rather obvious, since the 
density is outward decreasing, which is true provided that the equation of state is 
non-decreasing. 

Lemma 4.1. Consider an equation of state p = p{p) that is non- decreasing. Then 
U ^3/4 holds for all regular perfect fluid solutions associated with this equation of 
state. 

Proof. Consider the surface U = 3/4 in the state space and the flow of the system (jS]) 
through this surface. We obtain dCZ/dA] [7=3/4 = — 3T(1 + ct)(1 + 3a)V/6A ^ 0, since 
T > when dp /dp > 0. (One can say that the surface [7 = 3/4 acts as a "semi- 
permeable membrane" for the flow of the dynamical system, since solutions can pass 
through this surface in only one way.) Similarly, for a surface f/ = C/q > 3/4 we 
obtain dU / d\\u=Uo < 0- Since, for regular solutions, U 3/4 as A — > —00, it 
follows that [/ ^ 3/4 from below as A — > —00 and that [7^3/4 for all A. □ 
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4 Proof of the theorem 



We now proceed to give the proof of Theorem 13. li Note that (a variant of) the 
proof in the special case g = (the proper Buchdahl case) has been given already 
in [6]. Here we extend these techniques. 

Proof. Using that 



the function <I>, cf. (jl2p . can be regarded as a function on the state space. In 
the state space description, regular solutions emerge from the line of regular initial 
values {{U,V,i}) \ U = 3/4,1/ = 0, r2 G (0, 1)}. Evaluated in this limit, ^ becomes 
the function 

whose values range in [0,4). We obtain 



For each constant g G [0,4), the equation ^ = g defines a surface in the state space 
that divides the state space in two halves. We are interested in the subset 

$ > 5 , (26) 

since this represents the Buchdahl type inequality we want to prove. According 
to ()25p . the subset ^ > g contains the line of regular initial values up to a certain 
value of fi, which is given as the (smallest) value of ^1 such that a < (jg, where 
CTg = {3 — g + y/9 — 2g) / g. This means that, initially, i.e., for sufficiently small A (or, 
equivalently, for sufficiently small r), regular solutions with a < CTg satisfy ^ > g. 

Let us consider a regular solution with 

a < ^ - 1 (27) 

and the function $ along this solution, i.e., we track $(A) = $(C/(A), V{X), f^(A)) for 
A S (—00,00). Obviously, ^ > g holds initially (i.e., for sufficiently small A), since 
the r.h. side of ()27p is always less than ag. In Lemma [4. 2 1 we prove that the derivative 
of when computed on the surface <I> = g, is positive under the assumption ([2] 
when we also use that U ^ 3/4; i.e.. 



dA 



^ . (28) 

<i'=g 
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4 Proof of the theorem 



Therefore, since, first, ^{X) > g for sufficiently small A and, second, ^{X) cannot 
pass through the ^ = g barrier, it follows that 

^>g (29) 

for all A G M. Equivalently, ^ > g ior all r G [0,R) and $ ^ g for all r G [0,R], and 
the theorem is established. □ 

It merely remains to prove the following lemma. 

Lemma 4.2. The derivative of ^ on the surface ^ = g is positive for all U < 3/4 
if and only if 

a <^-l. (30) 
Proof. A straightforward computation based on the system ^ yields 



dX 



{3-AU){l-U + aU) 
= 3(1 -U)+aU - ' ^''^^ 



where ipg is given by 



^ 36{l-Uf -g (3(1- U) + aUf 

18{1 - U) {2{1 - U) + aU) - g {3{l - U) + aU)^ ' ^ ' 



Clearly, if ^7 < 3/4, 



d^ 
dA 



> if and only if ipg > . (32) 



Since the denominator of ipg is equal to the numerator of (pg (which we call ipg in 
the following) plus a positive term, namely 18(1 — U)aU, positivity of ipg already 
implies positivity of ipg. Therefore, 

^ >0 if <fg = 36{l-Uf -g{3{l-U)+aUf >0. (33) 
dX 

Positivity of (pg is in turn governed by a simple linear condition, 

ipg>0 ^ 6{l-U)-^{3{l-U)+aU)>0. (34) 

The derivative of the l.h. side w.r.t. U is given by —3(2 — y^) — a^/g, which is 
negative. Therefore, the inequality is satisfied for all U ^ 3/4 if and only if it is 
satisfied loi U = 3/4. Accordingly, 

^g>0 ^ ^(2-v^(l + a)) >0 (35) 
<^ 0- < — - 1 , (36) 
which proves the claim of the lemma. □ 
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5 Conclusions and outlook 



Remark. The lemma shows that the surface ^ = g acts hke a semi-permeable mem- 
brane in the state space (provided that U < 3/4 and a < — 1). Solutions that 
satisfy ^ > g cannot pass through the surface ^ = g and are therefore confined to 
the subset ^ > g of the state space. 

5 Conclusions and outlook 

In this paper we have proved a theorem that formulates a straightforward and simple 
generalization of the Buchdahl inequality for spherical models. Our considerations 
are based on the assumption of perfect fluid matter that satisfies a non-decreasing 
equation of state p{p); however, it is clear that the theorem can also be applied 
to collisionless matter (i.e., to the Vlasov-Einstein case), provided that one restricts 
oneself to isotropic configurations. We choose to not pursue the Vlasov-Einstein case 
further here; instead we will provide a comprehensive discussion of Vlasov-Einstein 
configurations (including some non-isotropic cases) and an associated generalization 
of Theorem 13.11 in another publication. 

In [7], assuming dominant energy, an inequality is derived, namely inequality (4.3) of 
that paper, which gives a bound on 2M/R by a number that is larger than 6/7. This 
inequality is claimed to be sharp; hence, in the light of Theorem 13. 1^ the solution 
with the maximal value of 2M/R given by [Tj (4.3)] must violate at least one of our 
assumptions. 

It is important to note that Theorem 13.11 makes no statement about whether the 
obtained inequality is sharp. (Trivially, it is sharp in the case g = 0, which is the 
original Buchdahl case, and in the case g = 4.) There might thus be room for 
improvement; whether this could be achieved based on the methods used in the 
present paper, by the techniques used in [7J, or a combination thereof, remains to 
be seen. 
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